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LONG-ITUDINAL STABILITY IN RELATION TO THE USE 
OF AN AUTOMATIC PILOT 
By Alexander Klemin, Perry A. Pepper, and Howard A. Wittner 

SUMMARY 

The effect of restraint in pitching introduced "by an_ 
automatic pilot upon the longitudinal stability of an air- 
plane has "been studied. Customary simplifying asTumptions 
have "been made in setting down the equations "of m otion an d 
the results of computations "based on the simplified equa- 
tions are presented to show the effect of an automatic pi- 
lot installed in an airplane of known dimensions and char- 
acteristics. The equations developed have "boon applied' "by 
making . calculations for a Clark "biplane and a Fairchild 22" 
monoplane . 

INTRODUCTION 



Pew theoretical studies of the automatic pilot in re- 
lation to the longitudinal stability of the airplane have 
been published to date. A number of interesting questions 
present themselves, of the following character: 



1. What is the effect on the longitudinal dynamic 
stability of an airplane of the gyroscopic stabilizing el- 
ement, introduced by the gyroscopic au tomatic - "pilot ? 

2. What modification of the longitudinal static s ta- 
bility is permissible when the gyroscopic stabilising ele- 
ment is introduced? — - 

3. How is the motion of the airplane in various typ es 
of gust affected by the introduction of a gyro scopic "sta- 
bilising element? 



4. Eow may the action of an automatic pilot be repre- 
sented mathematically? 

5. What is the effect of a lag in the action of tie 
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automatic pilot, so that the perfect gyroscopic stabilizing 
element is not available? 

6. What is the importance of the fact that, when an 
airplane runs into a vertical gust, some time passes before 
the tail surfaces are subjected to the influence of the 
gust? 

An attempt has been made in the present investigation 
to answer these questions by calculations applied rather 
fully to the Glark biplane treated in reference 1 and by 
preliminary calculations applied to the Fairchild 22 air- 
plane . 

The equations of motion have been presented in their 
simplest form, with a number of factors eliminated, to 
avoid intolerable complexity of calculation and to make 
the effect of the introduction of the automatic pilo,t more 
clearly evident. 



I. EQUATIONS OF MOTION WITH A GYROSCOPIC 
STABILIZING FACTOR 



The following equations are based on the system of 
axes used by the N.A.C.A. and the symbols have their usual 



significance, except for M u , M w , 
port these symbols are defined as 



w - du / 



m 



and Mg. In this re- 



etc. , but can be converted to standard form by dividing by 
k y 2 - 

If all derivatives are referred to unit mass of the 
airplane and M u is considered to be zero, the determi- 
nantal equation for horizontal motion becomes: 



I) - X u 

- z u 

0 



D - Z 



w 



w 



S 



= 0 



(I-D 
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If a motion in pitch away from the horizontal intro- 
duces a stabilizing factor Mg , the equation of angular 
motion "becomes 



3 



so that the de terminantal equation is now 



w 



- Z 



u 



D.~ Z 



W 



e 



U o D 



- Li 



w 



= 0 



Equation (i-l) may he written as 

A?\ 4 + B X 3 + OA 2 + D X + E = 0 



where 



B = - (X^ 2 + Z^E B 2 +.-M q ) 

G = X U M^ + X u Z w Kg 3 - X W Z U 2 B 2 + Z W M^ - M w u o 
D = X U M W U 0 - X u Z w li q + X w Z u M q 
E = Z u M w g - 
Equation (1-2) may he written as 

A l X 4 + B'X 3 + O'X 2 + D'X + E T =0 



where 



A i 
B t 
0 » 
D » 
E » 



= A 
= B 

= o - m 6 

= D + M 8 (X U + Z ff ) 

= 25 + MVu - X U Z W ) 



(1-2) 



(1-3) 



(1-4) 
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Derivatives for the Clark Biplane 

Because of- the completeness of the information avail- 
able and "because only the general effect of the introduc- 
tion of the factor Mg is to he investigated, the Clarlc 

biplane of reference 1 has been selected for experimental 
calculations.. The main characteristics and derivatives of 
this airplane are as follows: 

Wing area 464 sq. ft. 

Span (maximum) 41 f t . 

Stabilizer area . . 16.1 sq. ft. 

Elevator area 16.0 sq.. ft. 

Tail length ' 24,5 ft. 

Weight 1,805 lb. 

Radius of gyration in pitch .... 4.65 ft. 

Air speed 76,9 m.p.h. 

TT 0 '112.5 f. p. b. 

Wing setting , 0° 

X u -0.15 8 

X w 0.356 

Z u -0.57 

Z w -5.62 

M u , 0 

M w 3.2 - 

M q -192.0 

When the values of the stability-derivatives are in- 
serted in equation (1-3), there is obtained 
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so that 



A = 21.62 

B = 316.9204 
0 = 1,49 2.9 608 
D = 266.3290 
E = 58.7328 

A r = 21.62 

3' = 316.9204 

0 ' = 1,492.9608 - Mg 

B f = 266.3290 - 5.778 M Q 

E 1 = 5 8.73 28 ~ 1.0909 Me 

Tattle 1-1 gives the discriminant for various values 
of Mg , which, as will "be shown later, are attainable 
with aa appropriately geared, automatic pilot. 



TABLE 1-1 



Expressions for the Discriminant with Varying Values of Mg 

Me= 0 : 21.62X 4 + S16.9204X 3 + 1492.9608X 2 + 266.3290X + 58.7328 
" Me= -180 : 21.62X 4 + 316.S204X3 + 1672.9608X 3 + 1?06.3690X + 255.0912 
Mq= -360 ! 21.62X 4 + 316.S204X 3 + 1852.9608X 2 + 2346.4090X + 451.4496 
Me= -720 : 21.62X 4 + 316.9204X 3 + 2212.9608X 3 + 4426.4890X + 844.1664 
i: e =-1080 :21.62X 4 + 316.9204X 3 + 2572.9608X 2 4- 6506.5690X + 1236.8832 
M Q =-2160 :2l.62X 4 + 316.9204X 3 + 3652.9608> 3 - + 12746. 8090X + 2415.0336 



Variation of Routh's Discriminant with Ug 
?or the "biquadratic equation 



AX 4 + BX 3 + CX 2 + DX + E = 0 



Routh' s discriminant, whose magnitude is a rou^h measure 
of the stability, is 
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R 9 BCD - AD - B JE 

Eouth's discriminant has "been calculated for a number of 
values of M . The results are given in table 1-2 



TABLE 1-2 



M 9 


South' s discriminant 


0 


118. 6 X 10 s 


-180 


613 .0 


-360 


1 ,193 . 6 


-720 


2 , 596.0 


-1,080 


4,266. 1 


-2 ,160 


10 ,9 60 .0 



The results of this calculation are also indicated 
in figure 1. The introduction of Mq powerfully in- 
creases South' s discriminant. 



Variation in Hoots with Varying Values of Mq 

Tho roots for varying values of M@ have been calcu- 
lated "by G-raeffe's method, with subsequent improvement in 
accuracy obtained by E . B. Wilson's method, and the re- 
sults are given in table 1-3. Times tro damp to one-half 
amplitude have also been calculated in the usual manner as 
well as the periods, of oscillation in seconds. 

As might be expected on physical grounds, the intro- 
duction of Mg powerfully increases the dynamic stability. 

As Mg increases, tho period of tho short oscilla- 
tion decreases, and tho damping time lengthens somewhat. 
The long oscillation passes over into two aperiodic mo- 
tions, which are powerfully damped. .. - 



2J.A.C.A. technical Note Ho. 666 



7 



TABLE 1-3 





Root s 


. — 

Time to damp to one- 
half amplitude 

( 3 Q C . ") 

\ P u u « / 


Period of 
o scillation 
(see. ) 


0 


-7.2410 ± 3.74l4i 
•-.0834 ± .18191 


0.0957 
7. 8410 


1.6794 
34.5420 


-180 


-6.8847 ± 4.20851 
-0.3160 
- . 5 734. 


0.1007 
2.1935 
1„2088 


1.4930 


-3 60 


-6 .5130 ± 4.63201 
-0.2337 
-1.3989 


0.10 64 
2.9660 
. 4955 


1.3565 


-720 


-5.8625 ± 5.7536a 
-0.2127 
-2.7209 


0. 1182 
3 .2588 
. 25 47 


1.0920 


-1,0 80 


-5 .3480 ± 6.7270i 
-0 . 2062 
-3 .7564 


0.1296 
3 .3615 
.1845 


0 .9340 


-2,160 


-4.9235 ± 9.81911 
-0 . 2008 
-4. 6108 


0.1408 
3 . 45 19 
.1503 


0.6399 



The Order of Magnitude of Mq 

The order of magnitude of Mg may "be conveniently 

studied from the wind-tunnel tests of the H3.E.C.-2 report- 
ed in reference 2. This airplane is sufficiently similar 
to the Olark biplane to ser-ve as a "basis for calculation. 

For a "biplane model of 3.7-foot span, at an air speed 
of 40 foet per second \rith elevator neutral, the following 
pitching momonts were obtained: 

Pitching momont , M 
{in. —I'd. ) 

-0.05 7 
-.154 
-.226 



Angle 



of attack, 
(deg. ) 

0 
2 
4 



a 



so that 
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_9M „ _ o .0 53 



^degrees 

In t'lie elevator tests, with zero angle of attack, the 
following values were obtained: 

Elevator displacement , 8 e Pitching moment, M 

(deg.) ~ (in. -lb.) 

0 0.02 
5 -.49 

so that 

a_M m _ 0.094 

If it is assumed that 1° displacement in pitch would 
be followed by 3° displacement of the elevator, the rati^ 

3M /9*£ 3(0 .0 94) _ 5 3 
5B"- ' da 0.0 53 

so that a value of 315/30 = 6(3M/Bcx) is a practical possi- 
bility. 

Since 

so that M9 can be -calculated for various values of the 
ratio 

39 ' 3a 
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an / 
ae ' 


m 

da 


m 9 


0. 


5 


• -180 


i 




-360 


2 




-720 


3 




-1 ,080 


6 




-2,160 



II. SIMULTANEOUS VARIATIONS IN THE STATIC STABILITY 
AND THE GYROSCOPIC STABILIZING- FACTOR 
Calculations for the Clark Biplane 



By obvious manipulations, the constant for the "biquad- 
ratic 1 of the Clark biplane, at zero angle of attack, can 
"be written in the form 

A = 21.62 

B = 316.9204 

'• C = 1,132 - 112.5 M w - Me'. 

D = 209.2 - 17.78 1% - 5.778 Me 

E = -18.50 M w - 1.0 9 Me 

In order to obtain a reply to the question "What modi- 
fication of the longitudinal static stability is permissi- 
ble when the gyroscopic stabilizing element is introduced?" 
a number of calculations have been made in 'which both M w 
and Me were varied and the results are shown in table 
II-l. Values of these constants and for Routh'"s discrimi- 
nant are given for Mg = 0, -540, -1.080, and -2,160, with 
M w successively given the values 0, -0.8, -1.6, -2.4, -3.2 

-4.8, and -5.4. In table II-2 periods and times to damp 
to one-half amplitude have been calculated for represents^ 
tive combinations; B, is the period of the short oscilla- 
tion and is the time to damp to one-half amplitude of 

the short oscillation. Likewise T and T. v are the 

l a id 

times to damp to one-half amplitude of the two aperiodic 
motions replacing the short oscillation. for the long o s- 
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cillation, P 3 represents the period and T s the time 
to damp to one~half amplitude. Then T a and T a are the 
times to damp" to one-half amplitude of the two aperiodic 
motions replacing the short oscillation- Figure 2 gives 
a plot of Routh's discriminant for the same representa- 
tive cases. 

From these calculations and chart-g, it would appear 
that, once a large value of Mq has been introduced into 

the equations of motion, it would he a matter of indiffer- 
ence whether the static stability is double the normal 
value or reduced to zero. 



TABLE II-l 



Values of Routh's Discriminant for Changes in M w and Mq 





A 


& 


C 


D 


E 


BCD 


AD 8 


B 3 E 


South ' s 


















discriminant 


































M g = 0 








0 


21.62 


317 


1132 


209 


0 


75.1 XLO 6 


0.94X10 6 


0 


74.16XL0 6 


-.8 


21.62 


317 


1222 


223 


15 


87.3 


1.08 


1.50 *L0 a 


84.72 


-1.6 


21.62 


317 


1312 


237 


30 


98.6 


1.21 


3.02 


94.37 


-2.4 


21.62 


317 


1402 


252 


44 


112.0 


1.38 


4.43 


106.19 


-3.2 


21.62 


317 


14S2 


266 


59 


126.8 


1.53 


5.94 


118.33 


-4.8 


21.62 


317 


1672 


284 


89 


152.0 


1.75 


8.96 


141.29 


-6.4 


21.62 


317 


1852 


C313 


119 


184.0 


2.13 


11.99 " 


169.88 




M Q = -540 


0 


21.62 


317 


1672 


3329 


589 


17.65X10 8 


2.40tf.0 8 


5.92 KLO 7 


14.56XL0 8 


-.8 


21.62 


317 


1762 


3343 


604 


18.71 


2.42 


6.06 


15.68 


-1.6 


21.62 


317 


1852 


3357 


619 


19.70 


2.44 


6.22 


16.64 


-2.4 


21.62 


317 


1942 


3372 


633 


20.80 


2.47 


6.36 


17.68 


-3.2 


21.62 


317 


2032 


3386 


648 


21.85 


2.49 


6.51 


18.71 


-4.8 


21.62 


317 


2212 


3414 


678 


23.93 


2.54 


6.81 


20.71 


-6.4 


21.62 


317 


2392 


3442 


708 


2 6.10 


2.58 


7.11 


22.81 
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TABLE II-l -Continued 



M w 


A 


B 


C 


D 


E 


BCD 


AD S 


B S E 


Routh ' s 















discriminant 










M 6 = -1080 




0 


21.62 


317 


2212 


6450 


1178 


45.23 X10 


8.99 XlO 


11. 84X10 7 


35.06X10 8 


-v8 


21.62 


317 


2302 


6464 


1193 


47.17 


9 .03 


11.99 


36.94 


-1.6 


21.62 


317 


2392 


6478 


1208 


49.12 


9.07 


12.14 


38.84 


-2.4 


21.62 


317 


2482 


6493 


1222 


51.09 


9.11 


12,28 


40.65 


-3.2 


21.62 


ol / 


2572 


6507 


1237 


53.05 


9.15 


12.43 


42.66 


-4.8 


21.62 


Ol ( 


2752 


6535 


12 67 


57.01 


9.23 


12.73 


46.51 


-6.4 


21.62 


Ol r 


2932 


6564 


1297 


61.15 


9.32 


13.03 


50.53 






1 


&q = -2160 






0 


21.62 


317 


3292 


12669 


2350 


132.4x10 s 


34.7 X10 8 


236.1 X10 6 


95. 34X10 8 


-.8 


21.62 


317 


3382 


12683 


2365 


136.0 


34.8 


238.0 


98.82 


-1.6 


21.62 


317 


3472 


12697 


2380 


139.5 


34.9 


239.5 


102.20 


-2.4 


21.62 


317 


3562 


12712 


2394 


143.2 


35.0 


241.0 


105.79 


-3.2 


21.62 


317 


3652 


12726 


2409 


147.2 


35.2 


242.0 


109.60 


-4.8 


21.62 


317 


3832 


12754 


2439 


155.0 


35.3 


245.0 


117.25 


-6,4 


21.62 


317 


4012 


12783 


2469 


163.6 


35.4 


248.0 


125.72 



TABLE 1 1- 2. 



Periods P and Times to Damp to One-Half Amplitude, E 





M 

vr 






T 


l " " 


P 3 . 


T 3 a 


T 3D 


( sec . ) 


( sec . ) 


( sec . ) 


( sec . ) 


( sec . ) 


-2160 


0 


0. 775 




0.188 




0.097 


3 .45 


-6.4 


. 617 




.134 




.165 


3 . 45 


-1080 


0 

-6.4 


1. 210 
.792 


0 .176 
. 118 




0.104 
.246 


3.45 
3.45 


-540 


0 

-6.4 


2.040 
.938 


0 .124 
.106 




0 . 20.6 
.436 


3.45 
3 .48 










T 

la 


T 

ID 










0 




•0 .077 


126 




3 . 630 




0 


-3.2 


1. 679 




0 .0957 


34.542 


7.841 






-6.4 


1.130 




.095 




26.1 


7. 680 
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Calculations for the Fairchild 22 Airplane 

The calculations made on the Olark "biplane in hori- 
zontal flight have "been repeated for the Fairchild- 22 
high-vring monoplane , a more modern airplane of somewhat 
higher loading. In these" calculations a horizontal glide, 
power off, has been considered. 

The main characteristics and derivatives of the Fair- 
child 22 airplane are as follows: 

Wing area 171 sq-. ft. 

Span 32.83 ft. 

Stabilizer area 15.8 sq. ft. 

Elevator area 10.4 sq. ft. 

Tail length 14.69 ft. 

Weight 1,600 lt>. 

Radius of gyration in pitch ' 4.41 ft. 

Air speed 90.7 m.p.h. 

U Q 133. f-.p.s. 

Wing setting . -. 1° 

X u -0.0834 

X w 0.2290 

2 U : -0.4830 

-2.4850 

M u • ' * • • 0 

M w ....... -2.1300 

Mq. -64.0 

When these derivatives are inserted in equation (i-l), 
the expressions for the constants in equation (1-3) "become: 



n 



T " 
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A = 19.60 

B = 114.40 

0 = 170.33 - 133.0 - Mg 

D = 20.38 - 16.57 M w - 2.57 M 6 



E 



15 .77 M w - 0.318 Mg 



or 



A = 19 .60 
B = 51.40 - 

G = 6.23 - 2.57 M q - 133.00 M w - Mg 

- 0.32 U - 16.57 M w - 2i57 Mg 



D = 

23 = - 15.77 - 0.318 M e 



Effects of variations in M w and Me-- The effects 
of variations in M w and. Mg are illustrated in tables 
I 1-3 and I 1-4 and in the curves of figures 3, 4, and 5. 

The results of the calculations shown in tahle II-3 
and figure 3 have the same general character as the results 
of similar calculations for the Clark biplane. As Mg 
increases from zero to a value of —270, the long oscilla- 
tion is replaced "by two rapidly damped aperiodic notions. 



At the samo time, the p"eriod' of th.o short oscillation "be- 
comes more rapid and the time to damp to one-half ampli- 
tude lengthens somewhat. The Routh's discriminant increas- 
es to ro-latively enormous values as Mg incroases. 

There are very- clear lessons to "be derived from tahle 
II-4 and figures 4 and 5. When M w is given zero value, 
the period of the short' o scillation . "become s unduly long and 
the long o scillation .passe s over into two real roots, one 
of which is zero and the other is of small' negative value 
so that damping is small. But with M w = 0 and even a 
moderate value of Mg, the periods and the damping times 
"become, sati sfactory, and tho Routh's discriminant "becomes 
onornous in value. At other values' of M^.', it is the 
value of. Mq that has the predominating effect.^ 
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TABLE I 1-3 



Periods and Damping Times for E-22 Airplane; M w = - 2.13 

(original value) 



Me 


Short oscillation 


Long oscillation 


South ' s 
discriminant 


Period 
( sec . ) 


Times to 
damp to 
one-half 
amplitude 
( sec . ) 


Period 
( sec . ) 


Times to 
damp t o 
one-half 
ampl i tude 
(sec. ) 


0 


1.66 


0.241 


22.5 


12.52 


mm 


2.3 85 x 10 s 


-270 


1.26 


. 299 




.673 


3.54 


• 49 .44 


-540 


1.03 


.337 


mm 


.433 


4.40 


120 .60 


-1080 


. 817 


.372 


urn 


.348 


4.60 


335 .06 


-2160 


.&77 


.394 




.314 


4.93 


1049 .57 



It would appear that the value of M w , tha.fr- is, the 
magnitude of the static longitudinal stability, is imma- 
terial as long as a gyroscopic stabilizing elem-ent can be 
introduced. 



TABLE J 1-4 



Periods and Damping Times for P-22 Airplane 







Short oscillation 


Long oscillation 








Period 
( sec. ) 


Times to 
damp to 
one-half - " 
amplitude 
(sec.) 


Period 
(sec. ) 


Times to 
damp to - 
one-half 
amplitude 
(sec.) 


Bouth ' s 
discrimi riant 


0 


0 

-2.13 
-4.26 


14.00 
1.66 
1.18 


0.241 
.241 
.240 


22.50 
21.00 


12.52 
11.51 


5.31 


0.388 X 10 6 

2.385 

6.609 


-540 


0 

-2.13 
-4.26 


1.40 
1.03 
,863 


0.556 
.337 
.303 




0.214 
.433 
.627 


5.52 
4.40 
3.63 


73.25 
120.60 
169.97 


-1080 


0 

-2.13 
-4.26 


0.935 
.817 
.714 


0.465 
.372 
.336 




0.251 
.348 
.435 


5.52 
4.60 
4,3 6 


242.50 

335.06- 

428.62 


-2160 


0 

-2.13 
-4.26 


0.622 
.577 
.535 


0.441 
.39* — 
.362 




0.268 
.314 
.368 


5.52 
4.93 
4.79 


866.00 
1049.57 
1236.14 
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E ffects o f variations in a nd Ii q -with a constant 

value of 140.- .The interesting result that emerges fron the 

calculations shown in table II-5 and in figure 6 is that 
the South' s di scriminant and the roots of the equations 
remain satisfactory with a moderate value of 11^ even 

when llff = 0, provided that Mq = - 2,160. Only when "both 
M w and are zero, does instability appear and even 

then the unstable motion takes 9.85 seconds to double in 
amplitude . 

The preliminary conclusion is that, with the aid of a 
gyroscopic stabilizing element, both the static stability 
and the damping moment might be considerably reduced. 

Of course, if damping moment were considerably reduced 
by appreciably shortening the tail arm, there would be dif- 
ficulty in obtaining high values of Mg. 

nevertheless, a most interesting prospect of investi- 
gation and development appears to lie open, tending to de~~ 
signs in which the length of the fuselage would be consid- 
erably shortened. If full advantage were also taken of 
the possibility of decreasing the moment of inertia about 
the transverse axis, so that longitudinal restoring or 
controlling moments could be reduced, there would be an av- 
enue of approach to the tailless flying-wing type of air- 
craft . 

TABLE I 1-5 



Periods and Damping Times for I'-22 Airplane 









Short oscillation 


Long oscillation 






Mw 


Mq 


Period 
( sec. ) 


Times to 
damp to 
one-half 
amplitude 
(sec.) 


Period 
( sec. ) 


Times to 
damp to 
one-half 
amplitude 
(sec.) 


Routh's 
discriminant 


-2160 


-2.13 


-64 
-32 
0 


0.577 
.576 
.574 


0.394 
.759 
7.340 




0.314 4.93 
♦308 4.93 
.308 4.93 


10-i9.57 X 10 6 
549.11 
76.17 


-2160 


-1.06 


-64 
-32 

0 


0.601 
.596 
.595 


0.416 
.828 
46.000 




0*289 5.11 
.288 5.11 
.287 5.11 


958.8 
483.2 
38.21 


-2160 


0 


-64 
-32 

0 


0.622 
.617 
.616 


0.441 
.927 
a 9.85 




0.268 5.52 
.248 5.65 
.248 5.65 


866.00 
419.34 
.25 



a Real part was positive, 
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III. MOTION III A SHARP VERTICAL GUST WITH AITD 
WITHOUT THE INTRODUCTION 1 OP A GYROSCOPIC 
. STABILIZING FACTOR 
Mathematical Statement 



The motion following a sharp upward vertical gust, of 
magnitude w Q (where w Q is positive), is given by 

(D - Z u )u - X w w + g6 • = X w w 0 ^ 



+ gQ 

- Z u u + (D — Z w )w - U 0 D6 = Z w w 0 ^ (III-l) 

0 - M„w + (K B 2 D 3 - MgD ) 6 = K w w 0 



The equations assume that the line of action of the gust 
is along the Z axis. 



Let 



D - X 



u 



- X 



- z 



v. 



w 



D - Z v , 



g 



- U 0 D 



K-n S D 2 « ].I„B 



L w 



= PCD) 



come : 



B ~ - -q_- 

The symbolic solutions for the three variable o then be- 



u = w, 



V/ = w, 



e = w. 



x w 




- x w 






g 


z w 


D 


- K 








"w 




- w w 




D s 








F(D) 








D. - X^ 




x w 








2u 




z w 








0 




M w 


Kb 2 


D 2 








P(D) 








D - X u 




- x w 


x w' 






- z u 


D 


- Z w 


z w 






0 




- IV 










P 











(III-2) 
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The total acc eleratio n along the OZ axis is 



Dw - U 0 D6' 



(III-3) 



The act ual vel o cit y in • space along the OZ axis is 
given by 



w. 



- *u 



w 



'17 



1.L 



g 

- U c 3) 
K 3 3 D e - M q D 



F(D) 

The actual solutions may he obtained, from .the modified 
Heaviside identity 



F(D) ' F(0) A A F ' (X) • 



(II 1-5.) 



whose numerical evaluation has novr been dealt with in many 
publications. (See references 3 and 4=. ) 



H7hen the gyroscopic stabilizing factor Hq is intro- 
duced, the symbolic solutions assume the form: 



u = w c 









x w 




g 


z w 


D 
















M w 


K B a D 2 


- MqD ~ Mq 


- x u 










g 


- Z u 


D 




Z W 






0 








E B 2 D 2 





(III-6) 



with similar expressions for w and 9. 

Solutions of the Equations of Motion . 

M otion i n- w.- Solutions of the equations of motion 

are given in tables III-l, III-2, III-3, III-4, IlT-5/ana" 
III-6. 



TABLE III-l 
Solutions for »/w 0 



n 6 



-180 
-360 
-720 
-1080 
-2160 



Constant 
term 



-1.0000 



-1.0000 
-1.0000 
-l.QOQO 

■l.dodo 

•1.QD00 



Long oscillation 



-0.0884t 

-Q.0281e cos(0.18l9t - 1.4551} 



Slowly 

dsnrped 
exponential 



-0.0608e 



-0.3160t 



-.0148e 



-0,233 7t 



-.OllOe 



~0.2127t 



-.0075© 



-0.2062t 



•.0Q79e 



-0.3008t 



Rapidly 

damped 

exponential 



0.1055e 



,1339e 



-0.5734t 



-1.3969t 



.393 Be' 



,7209t 



.6659e 



-3.7564t 



.8395e 



-4. 6108t 



Short oscillation 



-7.2<il0t 

l,Q949e cos(3.74l4t - 0.4137) 



6 8R47t 

.9793e~ " cos(4.2085t - 0.2373) 

.8898e~ 6 * 5130t cos(4.6320t - 0.0808) 
,6496e~ 5 ' 8625t cos(5,7536t + 0.2599) 
,4578e" 5 * S480t coe(6.7270t + 0.4092) 
*1932e~ 4,9235t cos(9,8191t 4- 0.5076) 



CO 



« 

o 



1^ 
© 
o 

o 

p 

H 

tel 
o 

d- 

a 
o 



Oi 



TABLE I I 1-2 
Solutions for u/ir 0 



m 9 


Constant 
term 


Long oscillation 


Short oscillation 


0 




-0.0884t 

0.3810e cos(0.1619t - 1.4862) 


7 2410t 

-0.0356e " cos(3.7414t + 0.4258) 






Slcwly 
damped 
exponential 


Rapidly 

damped 

exponential 




-^160 




0.0768e--°- 2006t 


-0.0 7 «3e- 4 - 6108t 


0.0033e~ 4,9235t cos(9.8191t - 0.7477) 
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Table III-l gives the motion in w (i.e., along the 
OZ axis) for varying values of Mq. When the airplane en- 
counters the sharp vertical gust w 0 , directed vertically 
upward or, rather, negatively along the OZ axis, the sum 
of all the terms evidently represents the actual motion in 
space. If the constant term is eliminated, the rest of the 
solution represents the motion relative to the upwardly 
directed flow of air. As already indicated in section I, 
the introduction of Mq splits up the long oscillation 
into two exponential terms, one slowly damped and the other 
rapidly damped. Study of table tli-l indicates that, as 
Mq is increased progressively: 

1. The short oscillation "becomes less rapidly 
damped, out the amplitude of this oscillation dimin- 
ishes. 

2. The slowly damped exponential term "bec ome s 
progressively less damped, "but its amplitude becomes - 
more and more insignificant. 

3. The rapidly damped exponential increases in 
amplitude, "but its damping steadily increases^ 

If the airplane does entor a region where there is a 
steady vertical motion of air, then tho sooner it adjusts 
itself to this condition and the sooner the "residual" mo- 
tion disappears, the "better. The conclusion is therefore 
that increasing values of Hq progressively improve the 
motion in w. ' ~" ~ ~~ 

This conclusion is strengthened by a consideration of 
the curves of figure 7 where the residual motion Th w/w 0 
is plotted against time; for the extreme cases considorea, 
Mq = 0 and Mq = - 2,160. At t = 0, the residual veloc- 
ity w/w Q = 1 in both cases, of course, "but the negative 
values of w/w 0 that appear after a second or so are very 
much smaller in the case of the large value of Mq and the 
motion as a whole disappears quite rapidly. 

H otion in u .- Results of calculations of the extreme 
cases for motion in u are tabulated in table IH-2. As 
might be expected on physical grounds, no constant term is 
present in tho solutions of the equations of motion . As 
in tho case of motion in w, the long oscillation splits 
up into a slowly damped and a rapidly damped exponential 
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term when the large value of Mq is introduced. The short 
oscillation remains "but has a considerably shorter period 
and less damping, just as in the case of motion in w. 
The value of the constants being so small in all terms, it 
is clear that the amplitude of the motion in u is great- 
ly reduced "by the introduction of the gyroscopic stabiliz- 
ing- factor. 

M otio n i n. 6 .- The improvement in the pitching moment 
is still more marked, as indicated by the equations in 
table II 1-3 and the curves of figure 8. After a second or 
so, the motion in 8 practically vanishes. 



Solution of the Equations of Acceleration 
along the OZ Axis 

From the solutions of the equations of acceleration 
along the OZ axis, a marked improvement is evident with 
the resultant acceleration vanishing very rapidly, as is 
apparent from tables III-4, III-5, and III-6 and from 
figure 9. There is similarity in the effects produced on 
w, dw/dt , and U 0 q and hence on (dw/dt) -» TJ 0 q. In 
all these motions, the short" oscillation remains with more 
rapid period, a trifle slower— damping, but lesser ampli- 
tude. The long moderately damped oscillation splits up 
into two exponential terms: one is slowly damped but of 
very small amplitude; the other is of greater amplitude 
but very heavily damped. The reason for the improvement 
is apparent from the solutions for dw/dt and dG/dt. 
The values of dw/dt damp out very rapidly, and the val- 
ues of d8/dt are much smaller (as would be exp-ected 
from the small values of 8). Although the introduction 
of the gyroscopic stabilizer cannot possibly have an effect 
on the initial value of the acceleration, which is depend- 
ent only on the value of Z w w 0 , the virtual disappearance 

of the pitching motion serves to eliminate the centrifugal 
forces that subsequently form the larger component of the 
acceleration. 



■| ■ ■ 



TABLE III-3 



Solutions for G/w 0 


Mq 


Cone taut 
term 


Long oscillation 


Short oscillation 


0 

-8160 




-0,0884t 

-0.0224e cos(0.1819t - 0.4061) 


-7.2410t 

0.0049e cos(3.7414t - 1.0948) 

• 

-4.9235t 

0.00154e cos(9.8235t -0.0283) p 


Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 


-0.2008t 

O.OOOOle 


■ -4.6108t 
~0.00150e 




TABLE III-4 g 

/dw \, & 
Solutions for - U Q q )j* 0 P 


Me 


Constant 
term 


Long oscillation 


P 

Short oscillation H 

a 


0 

^2160 




-0.0884t 

-0.0477e , cos(0,18l9t - 1.4043) 


-7.2410t , . 3- 
-6.1235e cos (3. 741 4t -0.4117) a> 

o 

■ 

cn 
o» 

~4.9235t 

^1.08646 cos(9.8l91t +0.5058) 


'Slowly 
damped 

exponential ■ 


Rapidly 

damped 

exponential 


-0.2Q08t 

0,0012 e 


. ~4.6108t 
-4, 6734e 



TABLE III-5 
Solutions for ~/w 0 



Me 


Constant 
term 


Long oscillation 


Short oscillation 


0 


— 


_q 0884t 

~0.0057e * cos(0.1819t - 0.3535) 


-8.9239e" 


■7.2410t 

coa(3.7424t 


- 0.8907) 






Slowly- 
damped 
exponential 


Eapidly 

damped 

e aponential 












„ n«n * -0,2008t 

0.0016e 


„ nnnr> -4.6108t 
~3.8708e 


-2.1222 b" 


-4.9235t ,„ 

C0B(9.8191t 


- 0.5984) 








TABLE III- 
Solutions for 


-6 






Uq 


Constant 
term 


Long oscillation 


Short oscillation 


0 




— 0 0S84t 

0.0450e * coB(0.iei9t - 1.5133) 


-4.4888e" 


" 7 - 2410t sin(3.7414t 


- 0.0000) 






Slowly" 
damped 
exponeptial 


Eapidly 

damped 

exponential 




i ' 




-2160 




n MM -0.2006t 
-O.OOOSe 


n „ nf .„ -4.6108t 
0.7?63e 


-1.8563e' 


-4.9235t 

cos(9.8191t 


- 1.1317) 
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IT. MOTION IN A BELAYED GUST "WITH AND WITHOUT 
THE INTRODUCTION OF A GYROSCOPIC STABILIZING FACTOR 



Calculations for motion in a delayed gust of the. form 
w 0 (l - e-'k) have been carried out for only the two ex- 
treme cases, JJg s 0 and Mg = - 2,160. The solutions_ 

are presented In tables IV- 1 to IV-4 and in the curves of 
figures 10, 11, and 12. The results shoiT tendencies simi- 
lar to those present under the action of the sharp gust. 

1. With the introduction of a large value of Mg, 
the long cscillation again splits up into two exponentials, 
one slowly damped and one rapidly damped. _ 

2. The motions in w and u are very much better in 
the constrained airplane . 

3. The motion in 8 practically disappears on the 
introduction of Uq . 

4. The maximum acceleration along the OZ axis is 
slightly smaller when Uq = 0. Without the introduction 
of the gyroscopic stabilizing factor, the less constrained 
airplane (Mg =0) has a tendency to head into the result- 
ant wind and to relieve the acceleration. The difference 
is so slight, however, as to be negligible. Subsequently, 

the acceleration damps out very much, better in the con- 

strained airplane. 



TABLE IY-1 
Solutions for w/w 0 





• ■ — ■■ 
Constant 
term 


Exponential 
term 


Long oecillation 


Short oscillation 


0 


-1.0000 


1.14436"* 


-0.0302 e ~°* Q884t ain(0.1819t-0.0843) 


-0.1503e~ 7,241C>t cos(3.74l4t+0.l262) 








Slowly 
damped 
exponential 


Eapidly 

damped 

exponential 




-2160 


-1.0000 


1.2402e~ t 


ft nnoc, "0.2008t 
-0.0O99e 


-0.2325e- 4 ' a0Bt 


0.01B3e" 4 * 9235t sin(9.8191t+0.1263) 










TABLE IV-2 










Solutions for u/w 0 






Constant 
term 


Exponential 
term 


Long oscillation 


Short oscillation 


0 




0.04286*"* 


— Q 0fi84t 

-0.3686e cos(0 t 18l9t+1.4405) 


—7 24l0t 

O.Q047e cos(3.7414t-0,1470) 








Slowly 
dantpsd 
exponential 


Rapidly 

damped 

exponential 




-2160 




i 

-o.neie"" 1 


0.0948e- C ^° 08t 


0.02l8e" 4 * 6108t 


-O.QD04e"" 4 ' 92S5t cos(9.8191t+0.4638) 



no 



Sz! 

ft 
* 



o 
o 

t? 

o 
p 

i— 1 

o 

ct- 
CD 

o 



cr> 



i 



I 



TABLE IV-3 
Solutions for 9/w 0 



Me 


Constant 
term 


Exponential 
term 


Long OBdllation 


Short oscillation 


0 


0,00000 


0.00276*"* 


-0.0025e~°'° 884t cos(0.1819t-0.588O) 


~0.0007e" 7 ' 2410t cos(3.74l4t-0.4637) 








Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 




-2160 


0.0000 


-0,0003866"* 


n n^r ~0.2008t 
0,000015e 


0.000429e" 4 * 6108t 


-0 . 000429e" 4,9335t cos ( 9 . 819H+1 . 1593 ) 



fc4 



TABLE IV-4 









Solutions for (~ 






Constant 
term 


Exponential 
term 


Long oscillation 


Short oscillation 


0 




-0.83876"* 


0. 0414e*"° * 0B84t cos(0.1B19ttl.2513) 


-7 2410t 

0.8419e cos(3.7414t* 0.1271) 








Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 




-2160 




-1.283 ee"** 


O.OOlSa- 0 ' 20081 


1.2943e- 4 - 6108t 


-0 . 1029 e~ 4 ' 9235t cos (9. 819 lt-1 . 4470) 



HI 
o 
p 

O 



S3 
O 
ch 
CD 



CD 

cn 



IV) 



i 
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V. MOTION IN AN INFINITELY SHARP HEAD-ON GUST WITH 
AND WITHOUT THE INTRODUCTION OF A 
GYROSCOPIC STABILIZING FACTOR 



Calculations have "been made for the motion of the 
same airplane, the Clark "biplane, encountering an infinite- 
ly sharp head-on gust of strength u 0 . " at t = 0. 



(D - X u )u - X wW 



+ g6 



- Z + (D - Z^)w - u 0 d e 
0 



z u* 0 y (v-i) 



- M w w 



+ (K B s D a - H qL D)9 = 0 



The introduction of the gyroscopic element modifies 
the pitching-motion equation to 



0 - VLyVr + (K S 2 D 2 - M D - M 0 .) 9 = 0 



The solutions for the motions in w, u, and 0, and 
for the acceleration (dw/dt) - U 0 q are. given in terms 

of u Q in tables 7-1, V-2, 7-3, and 7-4 and , in the curves 
of figures 13, 14, 15, and 16. In figure 14 the residual 
motion in u is plotted, that is to say, with the con- 
stant term u 0 eliminated, "because the interest is mainly 

in knowing how quickly the airplane settles down to its 
final velocity U 0 - u 0 . 



TABLE V~l 
Solutions for w/uq 





Constant 
term 


Long oscillation 


Short 


oscillation 


0 




-0.0857e"°'° 8B4t cos(0.1819t + 0.4296) 


0.0779e- 7 - 2410t 


cos(3.?414t + 0.0737) 






Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 






-^160 




-o.weea- 0 - 2 * 8 * 


n 1Anc -4. 6108t 
0.1075e 


-O.WOle- 4 - 9235 * 


coe(9.8191t - 1.5150) 



S2J 

o 

* 



o 
o 

P 

H» 
O 

to 
H 

!4 
o 

cl- 
CD 

5*5 

O 



o 



TABLE 7-3 
Solutions for u/uq 





Constant 
term 


Long oscillation ■ 


Short oscillation 


0 


-1.00000 


1 . llS58e"°' ° 884t cos ( 0 . 1819 ttO . 4504) 


-0 . 002 54e~ 7 ,2410t cos ( 3 . 7414t +0 . 9045) 






Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 




-2160 


-1.00000 


l.OMSSe- 0 - 2006 ' 


-O.OlOOSe- 4 - 6108 * 


0.00017e~ 4 ' 9235t cos(9.8191t+0.3656) 



ro 
-a 



TABLE 7-3 



Solutions for 6/uq 



Me 


Constant 
term 


Long oscillation Short oscillation 


0 

-2160 




-0.0884t 

-O.Q0712e cos(0.1819t+1.5299) 


-7.2410t 

0.00035e cos(3.7414t-0.6081) 
O.OQ008e" * cos(9.8191t+1.0909) 


Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 


-0.2008t 

0.00016e 


n nnmn -4. 610 8 t 
-0.00020e 






I ABLE V- 
Solutione fo 


4 

r (aT " m)/^ 


Mq 


Constant 
term 


Long oscillation Short oscillation 


0 

-3160; 




-0.0884t 

-0.15361s cos(0.1819t+0.5165) 


-?.2410t 

-0.43684e cos(3,74t+0.0703) 
O.OSeeOe" 4 " 933 ^ cos(9.8191t-I.5l30) 


Slowly 
damped 
exponential 


Rapidly 

damped 

exponential 


« -0.2008t 
0.02478e 


rt ennAj -4.6108t 
~0.59844e 



■ « til 

jl ■ I I I IB I ■ I ■ I I im i i ■ ■ mm mn ■ ■ ■ ■ 
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As noted in the previous sections, with the introduc- 
tion of the gyroscopic stabilizing factor, the long os- 
cillation passes over into tiro rapidly damped aperiodic 
motions. The short oscillation has its period lengthened. 
The amplitudes "become very small in all cases. The im- 
provement in the motion due to the introduction of the 
gyroscopic stabilizing factor is greater than is the case 
with vertical gusts. This fact is readily explainable by 
simple physical considerations. In the case of the ver- 
tical gusts, in the first instants of the motion, it is a 
drawback that the gyroscopic stabilizing factor tends to 
maintain an invariable attitude relative to the horizontal; 
whereas, without the automatic pilot, the inherent stabil- 
ity tends to vary the angle of attack so as to reduce - the 
accelerating force along the normal axis. But in the case 
of the head-on gust, the action of the automatic pilot is 
beneficial from the very first instant, as quite elemen- 
tary considerations indicate. The incidence of the gust 
increases the lift and produces an upward motion that de- 
creases the angle of attack. In an ordinary"airplane tbXs 
action would be followed by a positive displacement in 
pitch. Such a positive displacement is vigorously resist- 
ed by the pilot, and hence the normal force and the ver- 
tical accelerations are kept down to a minimum. This re- 
sult is clearly indicated in figure 15, Indirectly "all 
other motions are improved. ■ ■■■ ■"■ - — - 



VI. MATHEMATICAL REPRESENTATION Off THE ACT IQET OE 
THE S FERRY AUTOMATIC PILOT 
Simplified Description of the Sperry Automatic Pilot 



Reference 5 gives a simplified, yet reliable descrip- 
tion of the Sperry Automatic Pilot. ' _ 

The functioning of this device can be followed by ref- 
erence to figure 17, which may be taken as applying "t~o~the 
elevator control. The rectangular box indicated in the up- 
per part of the diagram contains the gyroscope and air 

pick-off s. Mien the automatic pilot is operating the Tox 
is under suction, produced by the suction pump. Air, in-" 
dicated in figure 17 as entering through the bottom of the 
box, drive s the . gyro scope by impinging on buckets on*"the 
periphery of the gyroscope rotor. The axis YT « of the 
gyroscope may be considered as remaining truly vertical 
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■because o.f a mechanism incorporated in the instrument (hut 
not shown in the diagram), which is designed to correct 
any precessional movement of the gyroscope. 

As shown in figure 18, when the airplane is flying on 
the desired flight path, both air pick-off s are partly 
open and pressure is equal on "both sides of the air-relay 
diaphragm. When the airplane pitches, however, the air 
pick-off s rotate with it, while the gyroscope axis retains 
it-s vortical position in space. The relative angular move- 
ment "between the air pick-offs and the gyroscope causes 
one pick-off to he more or less completely uncovered while 
the other pick-off opening is blocked of f- a like amount. 
The result of pitching the nose of the airplane up is illus- 
trated diagramnat ically in figuro 17, where the. loft pick- 
off is shown fully open and the. right pick-off completely 
blocked off. 

Unequal restrictions are thus imposed on the flow of 
air entering the air relay through the air inlets on each 
side of the diaphragm and oxhausting into the gyroscope 
chamber, resulting in a displacement of the diaphragm ow- 
ing to the unequal pressure drop along the two air paths. 

The oil-relay valve, being mechanically connected to 
the air-relay diaphragm, receives the same displacement as 
the diaphragm, thereby admitting oil under pressure to the 
hydraulic servo-^cylinder.. The servo-cylinder piston, as a 
result, is displaced to the left. The piston rod, connect- 
ed to the servo-cylinder piston, protrudes tlrough both ends 
of the cylinder and, for the case under consideration, 
would be hooked into the elevator-control system in such a 
manner that a movement of the piston to the left would de- 
flect the elevator downward. 

A speed-control valve Z is used to regulate the 
speed of response of the piston in accordance with condi- 
tions of gustiness. A pressure regulator P prevents the 
speed of travel of the servo-cylinder piston from exceeding 
a certain value. 

Ah important element of the gyropilot is the follow- 
up system. The function of the follow-up is to insure 
that the control movement' will be proportional to the air- 
plane displacement. It is clear that, without the follow- 
up, any small displacement of the airplane would causo the 
el-ovator deflection to continue to increase until it 
reached the maximum value, unless the airplane in the mean- 
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tine had returned to the preassigned flight attitude. In 
practice, this characteristic would result in excessive 
control application, making the airplane overshoot when 
returning to the desired flight path. Thus any sli ght dis- 
turbance would cause a series of violent, rapidly diver- 
gent, oscillations. 

The follow-up system consists of a cable and pulley 
arrangement attached to the servo-cylinder piston rod so 
that the pulley S (fig. 17) on the gyro scope" case is ro- 
tated in proportion to the displacement of the servo— 
cylinder piston. At the other end of the shaft carrying 
the pulley S is a worm and wheel sector T, which thus 
rotates the air pick-off assembly in proportion to the 

displacement of the servo-cylinder piston. When the air- 

plane receives a large disturbance in pitch, considerable 
displacement of the servo-cylinder piston can occur before 
the air pick-offs are neutralized by means of tho follow- 
up system and further control displacement is prevented. 
For small disturbances, however, little control deflection 
will occur before the follow-up rotates the pick-offs back 
to neutral. 

Information Supplied by the Sperry G-yroscope Corporation 



In a communication received from H. Hugh Willis, Re- 
search Laboratory of the Sperry Gyroscope Corporation, 
the following information has been supplied whicTnTTs -- il- 
lustrated in figure 19. - 



The angular motion of the airplane that occurs 
before servo control starts varies from 10 to" "20 min- 
utes of a degree on commercial Gyropilots. Recently, 
we have found in the larger airplanes a dead spot of 
30 minutes of a degree provided better flying control 
and we are preloading the balanced oil valve accord- 
ingly. For bombing, mapping, etc., we adjust air 
gaps closer and obtain servo control for angular 
changes of 6 and 7 minutes of a degree. This causes 
continued hunting of the control at a frequency ap- 
proximately 250 cycles a minute. The amount of con- 
trol is insufficient to affect the attitude' of the 
airplane unless the airplane deviates more than 6 or 
7 minutes of a degree. 

. The speed of servo control varies proportionately 
with amount of deviation up to 2°. • The following 
curves show the relation of angular change of attitude_ 
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of the airplane with regard to the gyro and the re- 
sultant servo speed. 

The follow-up system limits the amount of con- 
trol per dogroo change of attitude of tho airplane 
as follows: 

For airplanes sensitive to control, the follow- 
up ratio is set to permit from 0.25° to 0.8° elevator 
movement per decree- change of attitude. Other air- .... 
planes require up t-o 2° control per degree change o f 
attitude. The average installation use s a 1:1 ratio. 

A linear motion of the servo-piston of 1 inch 
varies angular control motion from 4° to 15°. 

Further in actual flying the practice is to re- 
strict servo speed by the servo-speed control . valve s 
in the hydraulic system. In smooth air servo speeds 
are limited to from 1 inch per second to 1 inch per 
5 seconds effective on deviations greater than 2 . 
In turbulent air, speeds are increased up to maximum* 

Due to the fact that oil pressures are limited 
to permit tho human pilot overpowering the gyropilot 
in emergencies there are times in very violent air that 
control load exceeds servo power and no control oc- 
curs momentarily. 

Other pertinent information received from the manu- 
facturer is: 

1. The limit of the di-splacement of the eleva- 
tor by the piston is simply t"he usual stops on the 
control, cable. 

2. On return travel the curve of figure 19 also 
applies. 

3. TVhen the airplane has returned to the same 
attitude as the gyro, the servo-pi stx>n is, of course, 
at rest. 

Discussion of Mathematical Representation 

There is reason to believe that the curve of figure 
19 is derived from laboratory experiments, while the state- 
ments "For airplanes sensitive to control, the follow-up 
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ratio is set to permit from 0.25° to 0.8° elevator move- 
ment per degree change of attitude. Other airplane s~ re- 
quire up to 2° control per degree change of attitude." 
are derived from experience in flight. 

To give an exact mathematical representation of the 
action of the automatic pilot would offer insufferable 
difficulties, hut a qualitative discussion may he of in- 
terest and he sufficient for the purposes of this investi- 
gation. 

Let 8 e , he angular displacement of the elevator. 

t , time at which the air pick-off has "been 

sufficiently displaced for the automatic 
pilot to come into play. 

Then 



while 




6 Q dt 2 



e dt 



2 



Also, if E is the pick-off ratio, that is, the ratio 
"by which the displacement of the elevator, has to be multi- 
plied hjf in order to give the "backward" travel of the 
air pick-off and, if 6 is the lag angle before the x>ilot 

• - - = — " 

comes into play, then the servo-cylinder has a differential 
pressure acting on it as long as ■ - 



when 




the servo-cylinder ceases to operate and the elevator, to 
be displaced. 
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It is therefore clear that, if R (the pick-off ra- 
tio) is always substantially greater than 8, the action 
of the elevator will "be represented substantially "by the 
equation 

f 

a s e = e - e g 

Hence, it follows that the -stabilizing action of the 
automatic pilot can be represented with fair approximation 
by the term Mg(8 - B g ). 

From the information given, it would appear that 8 e 
is always a constant but," in reality, the automatic pilot 
provides for adjustment to suit varying atmospheric condi- 
tions, with piston acceleration greater under conditions 
of greater gustiness. 

It may also be surmised that the action represented by 
the equation R 8 e = 8 - 8 g is assisted by other flight 
conditions. When a vertical gust, say, first strikes the 
airplane, the angular acceleration is also at a maximum. 
After the angular lag has been overcome and before the el- 
evator has been displaced, the acceleration of the servo- 
cylinder should also be at a maximum. Inasmuch as the . 
angular acceleration of the airplane diminishes under the 
action of damping and elevator movements, the automatic 
pilot now encountering a greater hinge moment from the 
displaced elevator must also impose a lessened accelera- 
tion on the elevator. 

Further consideration indicates that, when the air- 
plane has received its maximum angular displacement so 
that q. = 0 , the elevator must also have its maximum dis- 
placement. Again, when the airplane "begins to return to .. 
its normal position, there must again be a lag in the ac- 
tion of the automatic pilot so that, its action must now be 
represented by the expression Mq(8 + 6 g ). 

711. MOTION IN A SHARP "VERTICAL GUST WITH 
LAG IN APPLICATION OP THE GYROSCOPIC STABILIZING FACTOR 



The motion in a sharp vertical gust , when there is a 
la S 6g in the action of the gyroscopic stabilizing ele- 
ment , can only he determined by a numloeT of calculations 
of a discontinuous character. 
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If, at time t = 0, the airplane is struck "by a sharp 
upward gust of constant velocity w Q , the equations of 
motion are those previously, given in section III, equation 
.(III-l): 



(D - X u )u - X^w + g e 

- Z u u + (D ~ Z w ) w - U 0 D8 



= ■ X W W 0 



= z 



W w O > 



(7II-1) 



- M w w 



+ (K B a D a - M a D)e = M^wo 



Thq motion represented "by these equations continues 
until a certain angular displacement 6g is reachod, which 

is equal to the angular lag of the automatic pilot. The 
value of 6 is determined by the size of the gap in the 

air pick-off (which is of the order of 7 1 as a minimum) 
and the combined or equivalent inertia of tho pilot. In 
these calculations 0g is assumed to "bo 0.5°, which prob- 
ably represents the maximum lag likely to bo used for a 
small airplane. ' " 

(D - X u )u - X^w + g 6 - x w w o 

- Z u u + (D - Z w )w - U 0 E 9 = Z w w 0 L(7II-2) 

0 - M w w 4- (K 3 s D a - Mq,D)e = M w w 0 

+ M 8 (6 - 6 g ) 

y 

The calculations for the solutions of equations (YII-2), 
in view of the calculations previously made in section III, 
were most conveniently made by splitting up the equations 
into two parts: 



(D - X..)u - X w w 



+ g e 



Z u u + (D - Z w )tt - U 0 V 6 
0 



- x w w o 



- v 



+ (K B 3 D S - M D 



K fl ) 6 = M w w 0 



and 
(D 



- X„)u - X. 



w 



w 



Z u u + (D - Z w ) w - TJ 0 D 9 



- M T ,w 



+ g G 

+ (K B 2 D2 ~ ic D - M Q )G 



_=.0 
= 0 



s 



> 
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In figure 20, curve A is the solution in $/vr 0 of 
equation (VII-4) and curve B is the solution in 8/w 0 ob- 
tained by adding the solutions erf equations (VII-3) and 
equations (VII-4) and so obtaining the solution af equa- 
tions (VII-2). In figure 8, the displacement 6/w 0 in 
radians per foot per second is plotted against time. As- 
suming' a value of 30 feet per second for vr Q and putting 

6 g = 0.5°, the value of 6/w 0 "becomes 0.00029. For the 

motion in which Uq 0, this value corresponds to a time 
tg = 0.15 second, approximately. 

After the lag has "boon overcome, the moti.on may he 
represented vory closoly "by equations (VII-2). In these 
equations the effect of the initial velocities u g , Wg, 

and q g at time tg = 0.15 second is neglected, because 

they are so small. It is evident that, owing to the lag, 
the gyroscopic stabilizing factor is indeed represented by 
the expression Mq(9 - 8 g ), where G g = 0.5/57.3 = 0.0087 

radian. Here, as in previous sections, Mq = — 2,160. 

Of course, the motion represented by curve S of figure 
20 or in equation (YII-2) does not continue indefinitely. 
The motion continues until a maximum displacement in pitch 
8max is obtained at which time q = dS/dt = 0 and t = 
t m = 0.25 second. (In reality, the actual time at this 

instant is t n + t g = 0.25 + 0.15 = 0.40 second.) 

At the instant when q = 0, the further displacement 
of the elevator ceases and, on the return displacement of 
the airplane, the action of the automatic pilot is repre-.. 
sented by the expression Mq(6 + 6g), owing to the lag. 
The reversal in lag in these calculations is taken as oc- 
curring instantaneously, which is not strictly accurate. 
To take accurate cognizance of the period of reversal 
would introduce overwhelming complexities, and no serious 
errors are introduced by this approximation because 0 g 

in practice will, for an airplane of" the size considered,, 
be considerably smaller than the 0.5° employed in thoso 
calculations. 

In order to allow for the reversal in lag", the solu- 
tion of the equations 
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(D - X u )u - X w w + g 8 - -=- 0 ^ ■ 

- Z u u + (D - Z^)w — TT 0 D 8 =0 1^(711^5) 

0 ' - M^w + (K B 2 D a - H q D-M )9 = 2M e e g 
is undertaken. . ^ 

The appropriate values for the solution of equations 
(VII-5) are plotted in curve 0 of figure 20. 

If curve C is displaced to the right to time t m = 
0.25 second and superimposed on curve 3, the resultant mo- 
tion shown in curve I) is obtained. This curve give~s~ ap- 
parently a constant value of 8 /w Q = 0.00025. But it must 
be remembered that, in the curves of figure 16, t = 0_ is, 
in reality, 0.15 second and that 6/tt 0 = 0 is, in reality, 
—0.00029. When appropriate changes are introduced for these 
initial conditions, the actual motion is that shown In fig— ~ 
ure 21, in which the fictitious constant value disappears. 

Table TII-1 contains a sumsfary of the pertinent solu- 
tions. A comparison of the curve of figure 9 with the 
curve of figure 21 (with Mg introduced) shows that iho 
lag of the automatic pilot does t change th.o character 
of the motion in 8 appreciably. 



TABLE TII-1 

For the equations 

(D - X u )u - X w w + g 8 = 0 _ 

- Z u u + (D - Z w ) w - U 0 D 8 =0 

0 - M w w + (K B a D s - MgD - Mg ) 8 = -Mg 8 g 

, -0.2008t -4.6108t 

8/w 0 = -0.000256 - 0.000036c + 0.000063o 

-4. 9 235 t ~ ■- 

+ 0.000269O cos(9.8191t - 0.3649) 
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(D ~ X u )u - X w w + g 9 

- Z u u + (D - Z w )w - U 0 D 9 

-0.2008t ~4.6l08t 
9/w 0 = 0.000512 + 0.000072e -0.000126e 

-4. 9 235 t 

- 0. 000538e coa(9.8191t - 0.3649) 

VIII. MOTION IN A VERTICAL G-UST WITH INTERVAL OF 
TIME EOR THE GUST TO REACH THE TAIL 
SURFACES TAKEN INTO ACCOUNT 



When an airplane passes into a region where there is a 
vertical gust, an interval of time elapses "before the ver- 
tical gust reaches the tail surfaces. The calculations of 
this section were made to ascertain the effect of - this lag. 



Elementary Theory of Gust Lag 

If the airplane enters a region in which a steady 
vertical gust of intensity w Q exists, there is an inter- 
val of time t = 1/U 0 (where I is the distance from the 
center of pressure of the wing to the center of pressure 
of the tail surfaces) during which tho tail surfaces are 
not affected "by the gust; the angle of attack of the tail 
surfaces relative to the wing is therefore le s aone d hy an 
amount w Q /U 0 . Thereby a stalling moment is introduced 
which, referred to unit mass, may he represented hy the 
equation 

where 



= 0 
= 0 

- M L - M Q ) 9 = 2M e 9 g 



r 
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P» 
m, 

3 t' 
I. 



density, 0.002378 slug/cu. ft. 
mass of airplane, 56 slugs. 

area of horizontal tail surfaces, 32.1 sq.. ft. 
distance "between centers of pressure. 



slope of lift of tail surfaces, equal to 0.0^3 
per degree. (This low value is due" to the_ 
low aspect. ratio of the tail surfaces, approx- 
imately 1.73.) 



w = 30 f .T5. s. 



so that 

M t = 127.0 ft .-lb. per unit mass. 
The interval of time t = l/U Q = 20/112.5 = 0.18 second, 



Motion with Gust Lag Taken into Account Without the 
Introduction of the Gyroscopic Stabilizing 

Element 

The motion with gust lag taken into account is the 
sum of the motions represented "by the two sets of equa- 
tions 



= ^rrW 



w»o 



(D - X u )u - X w w + g 6 

- Z u u + (D - Z w )tt - U 0 D 9 

0 - M^w + (K 3 S D 3 - l£ q D)'6 = M^w 0 



= Z w w o 



V (7III-2) 



y 



for which solutions have already been obtained in section 
III, and 

^ 



(D - X u )u - X w w 



+ g e 

- % u + (D - Z W - II DG 



u 



W' 



= 0 
= 0 



(YIII-3 ) 



- M^w 



y 



The motion in 8 without the vertical gust intervening 
derived from equations (VIII-3) is given 'by curve A of fig- 
ure 22 . 
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The motion in 8 due to the gust w 0 , obtained from 

the solution of equations (VIII-2) is given in curve B of 
figure 22. 

The resultant motion due to the gust w and to the 
stalling moment is indicated in curve G of figure 22. 

It is interesting to note that, although the vertical 
gust produces a negative displacement in G (as might be 
expected since the statically stable- aijplane tends to 
head down into the resultant wind), the stalling moment 
naturally produces a positive displacement in pitch. The 
initial disturbance of the airplane in pitch is there-fore 
ljess marked than when the vertical gust is taken as act- 
ing over, the entire airplane. 



Motion After the Effect of Gust Lag Has Vanished 

After the interval of 0.18 second, the vertical gust- 
embraces the entire airplane, and the ensuing motion must 
be calculated from the equations 



(D-X u )u-X^w + g 8 = Z w w 0 +Du o l 

-S u u +(D-Z w )w-U 0 D 6 = W^o f ~V 6 o ' 

0 -M w w + (K B 2 If -1^3) )9 = M w w D + K B 2 D S 8 Q '-M q D8 0 » 

+ K/Dq 0 » 

where - - 



y (tiii-4) 



u n 1 , is the velocity along the OX axis at t =» 0.18 
and u 0 f /^ 0 = 0.0330. 

w 0 ' , velocity along'the OZ axis at t = 0.18 and 
w o r /w 0 = - 0.5 720. 

8 0 ', angular displacement" in pitch at t =0.18 
and 6 Q »/ Wq = 0.00025. 

q ' , pitching angular velocity at t = 0*18 and 
q 0 '/w 0 = 0.0057. 

The. calculation is continued by introducing u Q 1 , w 0 ' , 8 Q 1 , 

and q 1 as initial values, 
o 



r 



V 
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In order to solve equations (YIII-4), it is conven- 
ient to separate them into the equations 



(3D - X u )u ~ X^w + g 6 . 

- Z u u + (D - Z w )vr - U 0 D 8 

0 - M^w + (K B a D 2 ~ Mq_D)e = M^w 



= X w*o 



for which solutions are already available from section III 
and the following equations, which allow for the initial 
conditions: 



= Du Q 1 



(2-X u )u-X TT w +g e 
-Z u u + (D-Z W ) W -XJ 0 3D G 

o ~m w w + (e b 2 if -m^d ) e = k b 3 d 3 e 0 e 0 i 



\ 



■ ^o' 



- u o D e o 1 



> 



(VIII-5 ) 



The solution in G from equations (YIII-5) is 

-7.241t 



6 / w o = " 0.00 25 6e 



cos(3. 7414t ~ 0. 8185 ) 



~0.0884t 

+ 0.00234e cos(0.18l9t-0.544l) (VIII-6) 

The complete motion, including that indicated in curve 
0 of figure 22 and the subsequent motion obtained by the " 
solution of equation (YIII-4), is shown in figure 23 . : _ The 
curve has been plotted only up to the time t = 3 . o' jsecohds 
from the initial instant when the vertical gust. first 
struck the airplane. 

It is interesting to compare thenotion in G of fig- 
ure 19 with that indicated in figure 8. When the "upward 
vertical gust is taken as striking the entire airplane in- 
stantaneously, the maximum angular disturbance is a "nega^ 
tive one. When the delay in reaching the tail surfaces is 
taken into account, the maximum angular disturbance "is ~a~ 
positive one, although it is much smaller numerically than 
when the delay is neglected. 



From this result, it is inferred that: 
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1. The vertical acceleration is certain to 'be great- 
er when the gust lag is taken into account, because the 
angle of attack of the main wings is actually increased by 
the action of the gust. " .7 

2. It has been sufficiently demonstrated "by previ- 
ous sections that the action of the gyroscopic stabilizing 
element is always to reduce the angular disturbance. Be- 
cause there is a greater angular disturbance with this 
type of gust lag, the automatic stabilizer should bo even 
more valuable under such circumstances. 



CONCLUSIONS 



Owing to the complexity of the subject, the present 
investigation can only be regarded as a preliminary one. 
Also the investigation suffers from the fact that calcula- 
tions have been restricted substantially to one airplane. 
Nevertheless, a method of approach has been established, 
and certain tentative conclusions may be sot down. 

1. The introduction of a gyroscopic stabilising ele- 
ment is, in general, a most powerful method of increasing 
the stability of an airplane. 

2. The introduction of a gyroscopic stabilizing ele- 
ment has a slightly detrimental effect on the short oscil- 
lation, whose period becomes more rapid and whose damping 
time increases. In an ultrasimplif ied conception of the 
short oscillation, velocity along the X axis may be con- 
sidered to remain constant. With this assumption the air- 
plane may bo imagined as running into a region in which a 
vertical upward current exists. Immediately the airplane 
tends to rise owing to the increased Qngle of attack but 
tBnds at the game time to nose down into the relative wind. 
Also, as the airplane tends to rise, the relative angle of 
attack again diminishes. Therefore tho" initial excess of 
lift does not exist very long, and the motion is rapidly 
damped, When, however, the gyroscopic stabilizing element 
is introduced, it tends to oppose the nosing down into the 
relative wind, so that the initial excess lift is somewhat 
greater than with' the naturally stable airplane. This ex- 
planation gives in rough and simplified fashion the detri- 
mental effect on the short oscillation. At the same time 
every calculation presented in this report indicates that 
the effect of the gyroscopic stabilizing -e-lement on the 
short oscillation is not serious. 
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3. The long oscillation may "be very roughly consid- 
ered in the. light of Lancbester's treatment of the phugoid. 
The airplane of natural inherent stability is supposed to 
nose down immediately into the relative wind and hence to 
have no excess lift. It therefore moves downward acquir- 
ing excess velocity and excess lift, until the excess lift 
gives it an upward velocity in excess of the upward veloc- 
ity of the gust. But this motion takes an appreciable' 
time; and the oscillation therefore has a slow period and 
slow damping. When, however, the gyroscopic stabilizing 
element is introduced, the travel along a downwardly in- 
clined path is soon checked, the attitude to the horizon- 
tal remains almost unchanged, and the excess lift soon dis- 
appears; . the airplane assumes the steady notion of the up- 
ward current of air, with rapid disappearance of the dis- 
turbance. This discussion is oversimplified hut explains 
the replacement of the long oscillation "by two well-damped 
aperiodic motions, as always indicated in the calculations 
of this report. 

4. When the coefficients 0, D , B of the "biquadratic 
equation are studied in detail, it is seen that the effect 
of the introduction of Mq overshadows the effect of It^. 
It may almost he. said that it is immaterial what the value 
of M w may he, provided that a realizable value of Mg 

ic introduced into the airplane. 

5. When, as for the Fairchild 22 airplane, the val- 
ues of 11 are varied, it is seen that Mq again over- 
shadows the effect of and, with considerably reduced 
values of M^, it is again possible to secure satisfac- 
tory values of periods and damping times by the use of "fig. 

6. Even when M w and are simultaneously re- 
duced, Mq proves a cure-all. _ 

7. It follows therefore that,, if an adequate gyro- 
scopic stabilizing element is introduced, the airplane de- 
signer need have much less worry regarding longitudibal 
static stability or damping. Interesting departures from 
classical airplane design are thus in prospect with the 
use of the automatic pilot. Of course, conservative air- 
plane designers will not at once abandon their efforts to 
secure inherent longitudinal stability; nevertheless, 
their thoughts ought well be turned in this direction 
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8. When the actual .mo' t ion of' tho airplano is studied 
with the introduction of tho gyroscopic stabilising ole- 
nont , it is seen that tho improvement in notion night "bo 
prodicted from tho .study of the biquadratic equation and 
its roots'. Such an inprovenent is found to occur in every 
type of disturbance studied? sudden vertical gust, delayed 
vertical gust, sudden head-on gust. It is in the decreased 
amplitude of the pitching motion that the improvement is 
most marked; again, when encountering the head-on gust, 

the improvement in the motion is even more marked than in 
encountering the vertical gusts. 

9. The mathematical representation of the action of. 
the automatic pilot is far from complete f Novcrthe.loss, it 
should "bo a close approximation to actual conditions. 

When the lag of tho automatic pilot is taken into account, 
with a magnitude greater than that likely to occur in 
practice for an airplano of tho size considered in this 
report, there is a lesser improvement than when the "ideal" 
gyroscopic stabilizing eloment is introduced. The charac- 
ter of tho improvemont in notion ronains . substantially 
tho same and, provided that the designers of the automatic 
pilot keop lag within close linits and prevent hunting 
within tho unit itself, it is fair to think of the auto- 
matic pilot as boing. close to tho action of the "ideal" 
gravitational stabilizer. 

10. Tho effect of delay of a vertical gust in reach- 
ing the tail- surfaco is shown to bo relatively unimpor- 
tant and not likely to affect the general conclusions. 

11. There is ample room for further investigation of 
the subject, which the tremendous potentialities of the 
autonatic pilot would fully justify. 



Daniel Guggenhein School of Aeronautics, 
New York University, 

New York, N.Y., June 1938. 
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Figure 7.- Residual w/w Q against time, t for constant-velocity 
gust w Q . 
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Figure 8.- Pitch angle 9 against time, t for constant- velocity gust 
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Figure 9.- dw „ against time, t for constant-velocity gust w a 
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